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ABSTRACT 

We consider particle oscillations and their damping in second-quantized form. We find 
that the damping or "decoherence" may be described by a Boltzmann-like collision 
integral with "non-abelian blocking factors" (fermions). Earlier results are generalized 
in that the momentum degrees of freedom are included and that the mixing equations 
become intrinsically non-linear at high densities. 
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^ Based, in part, on work to be submitted by G.S. as a doctoral thesis to the 
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1. Introduction. — A series of interesting problems, ranging from the question of 
"decoherence" , such as for the two states of an optical isomer spin relaxation in 
condensed matter^, or the damping of particle oscillations^, pose essentially the same 
question in a variety of contexts^. For neutrino oscillations in particular, questions 
of statistics can arise, as in the neutrino degeneracy during the collapsed phase of 
a supernova. It then becomes of interest to examine a second-quantized approach 
to such problems, where giving the correct commutation relations to the fields takes 
care of "Pauli blocking" and similar effects. In addition it might be hoped that 
the use of field theoretic language may add some insight into the interesting quantum 
mechanical issues behind our practical applications. In attempting such a formulation 
we have found a simple and transparent kinetic equation, which may be thought of 
as a kind of "non-abelian Boltzmann equation" . Although closely related or partial 
versions of this equation seem to exist in various connections^, we think it useful to 
present it here in a general and simple form. 

An interesting aspect of the equation is that it presents kinetic and "interfer- 
ence" aspects simultaneously. Spin-like problems of the type under discussion depend 
very much on subtle phase effects. On the other hand the Boltzmann type of kinetic 
equation seems to be entirely "incoherent" and classical in spirit. Nevertheless, we 
shall see how a marriage of both aspects arises naturally. This also may be of help 
in understanding many questions including how the classical world arises from the 
quantum mechanical one. 

2. Matrix of densities. — In a previous paper we gave a field-theoretic approach 
to mixing problems dealing directly with the fields. However, incoherence or damp- 
ing was not considered. In the single-particle framework used previously^ to study 
damping, the object of study (for two-state systems) was the quantum mechanical 
2x2 density matrix for the two states. The two-state system was viewed as imbed- 
ded in a much larger system and the density matrix was found by considering the 
total wavefunction of the whole system and tracing over the unobserved variables of 
the large system — usually thought of as "the medium". If the density matrix was 
expressed in terms of Pauli matrices r and a polarization vector P as ^(1 -|- P • r) 
then P was found to obey a Bloch-like equation: 

P = VxP-DPt. (1) 

This gives a coherent rotation of the polarization through V, and through the damp- 
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ing parameter D, a loss of coherence or damping manifested as a shortening of P. 
The subscript T indicates those components non-diagonal ( "transverse" ) with respect 
to interactions with the medium. Also, a formula for D was obtained in ^'-matrix 
language, resembling a generalized optical theorem, and giving D in terms of the 
scattering amplitudes on the background. In measurement theory language, D may 
be thought of as giving the rate of "reduction" of the wavefunction of the two-state 
system^. 

What is to play the role of the quantum mechanical 2x2 density matrix in 
the field-theoretic formulation? In second- quantized terms the above procedure of 
tracing over the unobserved variables may be accomplished in a sense by applying 
an annihilation operator to the state to "remove" the coordinate and then finding 
the overlap with another such operation. Thus we propose to study the "matrix of 
densities" 



instead of the density matrix. The indices refer to the various components of the field 
in question and | ) to our original state. For a two-state system, such as neutrino 
fields with two flavors, ip will be a two-component spinor in flavor space^, and we 
will have four bilinears to deal with. Diagonal elements of this c-number matrix are 
the ordinary densities of a particle type, while the off-diagonal elements contain more 
subtle phase information relating to the coherence of the system. The matrix (2) is 
Hermitian and positive definite. If we assume interactions that conserve the total 
number of particles, as we shall in this paper, then the trace of its spatial integral 
is conserved. Hermiticity, positivity and conservation of the trace are the properties 
characterizing a density matrix and so it is justified to speak of the integrated p as a 
density matrix^. 

The object (2) is of course much more complicated than the density matrix 
of the single particle formulation since the arguments of the fields can refer to any, 
and in general different, points of space; similarly in momentum space they can refer 
to different momenta. However, since our major interest is the mixing question, we 
shall make the greatly simplifying assumption of ignoring spatial inhomogeneities. 
Our state, | ), will represent a situation which is spatially uniform. Thus in position 
space p can only depend on the difference of the arguments of the fields, while in 
momentum space we have ( |7/^j(p')-?/'i(p)| ) = (27r)^5^(p' — p)pij(p), and so we can 
confine ourselves to the study of the matrix pp. (When matrix indices are suppressed 
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we write p as a subscript). 

In this paper we shall limit ourselves to the non-relativistic domain in the sense 
that we do not directly consider the role of creation and annihilation of particles. 
Thus our fields are not those of fully relativistic quantum field theory; the operator 
■0(x) contains only annihilation operators for particles and not the creation operators 
for antiparticles. 

3. Interactions. — Turning to the time development of p, we treat the fields 
as Heisenberg operators and assume that our constant state | ) is produced 
by letting the creation operators for the -0 particles operate on the state vector for 
the medium at i = 0. This corresponds to a simple product wavefunction with no 
correlations between the ^-particles and the medium. Wc follow the approach of 
finding the equations among the py by turning on the interactions for a short time 
at t = and then assuming these equations to be valid at all times. 

Our starting point is the equation of motion for the operator ^/'j^/'i, the Hamil- 
tonian consisting of a "free" and an interaction term, H = i/fj-cc + -f^int- -f^frcc is taken 
to be bilinear in the fields and diagonal in momentum. It may describe a kinetic 
energy plus some internal interaction of the fields among themselves. An example 
is the free Hamiltonian for neutrinos with a mass matrix, which leads to a rotation 
in the internal fiavor space. If refer to different spin components of a system 
with a magnetic moment, i^free may refer to an applied magnetic field, also leading 
to a rotation among the components. Or for the molecular isomers i^free gives the 
tunneling between left and right configurations in addition to the kinetic energy. 

iffree be written for each momentum in terms of a matrix Q, and the column 
vector ijj as ■0pQp'0p. Commuting with ■0j(p)'0i(p) produces bilinears of the same 
momentum. The set of densities "rotate" among each other. This is of course just 
the precession of the "spin" , ordinary or flavor variety. 

i^int gives the interaction with the background medium which we take to be 
a local point-like interaction between the ip and the background, represented by an 
operator B: 



G is a coupling constant matrix acting on the column vector of the ■^'s, and B would 
be typically, as in the Fermi interaction, bilinear in the background fields. When G 
is diagonalized it gives a set of states which scatter only into themselves. The T for 
"transverse" in Eq. (1) then refers to those bilinears which are off-diagonal in such 
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states. The choice of Eq. (3) is appropriate to the Fermi interaction and many other 
problems where the particles under study undergo a local elastic scattering as their 
basic interaction. We shall take G Hermitian, corresponding to lowest order, but one 
can imagine relaxing this restriction to represent an underlying process in higher than 
Born approximation, as in the ^'-matrix approach^'^. 

The commutator of i^int with ■0j(p)'0i(p) again yields a bihnear. Finally taking 
the expectation value of the resulting equation of motion yields [notation: dp = 
d3p/(27r)3]: 

ipy(i,p) = [Q(p),p(i,p)]y+ j dp'(|s(t,p-p')V'](i,p)G'iA:V'fc(*,p')|)-h.c. (4) 

This is an exact expression in terms of Heisenberg operators. It resembles but is not 
yet in the desired form of a set of equations among the py, hke Eq. (1), because B 
cannot be factorized out of the brackets. The original state | ) is given by free particle 
operators, ^/^^ at t = 0, while the operator expressions in the brackets are non-trivial 
in terms of such Schrodinger operators. 

We can, however, introduce perturbation theory where we expand the Heisen- 
berg operators il){t) and B{t) in terms of non-interacting (interaction representation) 
operators ipQ{t) and B{){t) corresponding to G = 0. A first approximation is to set 
il){t) = ipo{t) and B(t) = Bo{t). Now with this approximation the bracket in Eq. (4) 
may be factorized since by assumption the original state | ) has no correlations be- 
tween the 7/;-particles and the background. Only bilinears with equal momenta survive 
by virtue of the spatial uniformity represented by | ) and we obtain a closed equation 
among the p: 

iPp = Pp] + Pmed [G, Pp] . (5) 

Here, Pmed = ( 1-^(0) 0)1 ) is the medium density, assumed to be stationary and taken 
at momentum transfer zero. This relation will be recognized as adding an "index of 
refraction" or medium- induced energy contribution to the free mixing given by i^free- 
Note this index of refraction effect actually refers to differences in the index (as occur 
for example as in the MSW effect) and does not appear for a one-component field. 

In the two-level [or more generally SU(2)] case where fl, p and G can be 
expanded in terms of Pauli matrices [more generally SU(2) generators], e.g. pp — 
|np(l -|- Pp • r), with coefficients Vfree, P and G respectively, Eq. (5) becomes 
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Pp = Vp X Pp with Vp = Vfree(p) + Vmed where Vmed = Gpmed- This is the 
first term of Eq. (1) and agrees with the old result^, to first order in G. 

4. Incoherence and damping. — Now, Eq. (5) still contains no sign of any loss 
of coherence, i.e. any shrinking of P. This is understandable, since to this order of 
approximation there is no excitation of the medium. Only forward scattering has 
been taken into account, as is indicated by taking B at momentum transfer zero. 
If the medium does not change its state, no incoherence can be expected since the 
overall wavefunction, assumed to be a simple product at t — 0, has remained so. As 
is familiar in perturbation theory, the energy but not the wavefunction changes in 
first order. Thus to see the first evidence of damping or incoherence we must go to 
the next order. 

We stress that spatial homogeneity excludes sources of damping or incoherence 
due to spatially varying conditions. For example, if the magnetic field applied to an 
ensemble of spins varies in space, or if there arc inhomogeneities in the medium in 
which our neutrinos propagate^^, there will be an averaging or loss of coherence in 
the properties of the ensemble. In general a loss of coherence can always result from 
averaging over some parameters characterizing the system. However, here we only 
wish to address the true quantum mechanical loss of coherence. 

We go to the next order in G by expanding the Heisenberg operators to first 

order: 



^{t) = ^o(t) + I f dt' [ifint(t - t'), 


B{t) = BQ{t) + I f dt' [H;^t{t - t'), Bo{t) 





(6) 



To evaluate the new terms in Eq. (6) we must now make our assumptions on the 
nature of the system explicit. We assume that the interactions described by i?int 
can be taken as individual, isolated collisions where the particles represented by the 
ip go from free states to free states as in ordinary scattering theory. Thus the time 
for a single coUision is taken to be short compared to the time between coUisions; 
furthermore the very small energies represented by the mixing effects in Eq. (5) are 
neglected in calculating the effects of the collisions. Under these assumptions we can 
carry out what amounts to a Golden Rule calculation on the level of the fields and 
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find 



iP{t, p) = Vo(^, P) + ie~'^^^ J dp'duj5{Ep, -Ep + uj)Bq{uj, p - pOG'V'oIO, p'), 
B{t, p - p') = So(i, p - p') + « j d({d(\dw5{E^ -E^-u)x 



X 



i?o(a;, q' - q), i?o(t, P - P ) ^J(0, qOGVo(0, q) 



The 5 function arises from letting the time integrals go to infinity. We have dropped 
an associated principle-part integral since this has to do with a second order energy 
shift in the medium, which we already have to lowest order in Eq. (5). 

We may now substitute Eq. (7) into the r.h.s. of Eq. (4). By construction the 
state I ) has the property that only bihnear correlations of the type ( iV'Vl ) ^-^e 
present, so in carrying out the contractions in the brackets only p's will arise. It must 
be recognized, however, that an assumption is involved here, similar to that involved 
in the classical Boltzmann equation. While we construct the state at t = such that 
no other correlations are present, the philosophy of the method that the equation 
found by turning on the interactions for a short time at t = can be used for all 
times tacitly assumes that no other kinds of correlations can build up significantly. In 
this sense our derivation is Boltzmann-like and, for example, we exclude the possibility 
that our system might have important pair correlations. 

Carrying out, then, the contractions with only the survival of ordinary density 
bilinears, and evaluating all quantities at t = 0, we arrive at our main result: 

Pp = -i [f^p, Pp] - ipraed [G, Pp] 

+ / dp' [S{v' - p) (^-Pv)Gp,'G-,h.c. _ _ g(l-Pp0^gPP + h.c. j 

(8) 

where -S'(A), a function of the four- vector A, is the usual dynamical structure function, 
defined by jl^ dt e*^o*( |So(t, A)So(0, -A)| ). This is the desired generalization 
of Eq. (1) under Boltzmann-like conditions. It resembles the ordinary Boltzmann 
equation with the usual drift term, absent because of spatial uniformity, replaced by 

the commutators for coherent rotations, which give the free motion in "spin space" . 
The principle new features of Eq. (8) relative to Eq. (1) is the inclusion of the mo- 
mentum degrees of freedom and the nonlinearity in p introduced by statistics. The 
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main difference with the ordinary Boltzmann equation is the matrix structure of the 
"coUision term". The "loss" term contains such effects as the creation of polariza- 
tion by absorption while the "gain" term gives effects like the transfer of polarization 
through scattering. 

We have taken the original state | ) to be pure. If it is itself characterized by a 
density matrix then our various bracketed expressions can be interpreted as thermal 
or other averages. Similarly, if there are several species of background particles, 
uncorrelated with one another, their contributions can be summed separately in the 
collision integral. Eq. (3) is the simplest possible interaction. However, as long as 
the essential feature that the interactions of the different ijj with the background 
are identical up to a coupling constant is retained, we anticipate no change from 
more general interactions except for a generalization of the structure factor S. For 
example, when the spin structure of the background is accounted for S{p — p') will 
be replaced^^ by the probability W{p,p'), which is no longer simply a function of the 
difference {p — p'). It may be verified that this has no effect on our arguments. 

By integrating over p we may also study the total density matrix, with its 
associated total polarization: 



We see here from the commutator structure that the total number of particles, given 
by Tr(p), is explicitly conserved. Observe that in the special case where all types scat- 
ter with the same amplitude on the medium, i.e. G proportional to the unit matrix, 
the collision term in Eq. (9) vanishes identically. There is then no dynamical damp- 
ing, however there remain "de-phasing" effects due to the momentum dependence in 
the precession of the different modes. 

5. Limiting cases and applicaUons: Ordinary collision integral. — The simplest 
limiting case obtains when we have a one-component ip (or equivalently many compo- 
nents but no mixing terms, or mixing terms which can be simultaneously diagonalized 
with G) and the background is simply free particles, i.e. -B(x) = 0t(^x)0(x), 0(x) a 
free field. Then: 



where nmed(Q) ^-^e the occupation numbers in the medium (± for bosons or fermions). 



P = -«Pmed [G, p]-i j [f^p, Pp] dp + i / S{p- p') [G, (ppG(l - PpO - h.c.)] dpdp' . 



(9) 




(10) 
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In the low density limit where the blocking factor can be neglected, with a single 
medium mode Qmed) ^(q) = i'^'^)'^ Pmed^^ ~ Qmed)) S becomes simply the energy 
S function times 27rpnied- Because J dp'S is the scattering rate, Eq. (8) has the 
appearance of standard rate-of-gain and -loss terms. With only a single mode p oi t/j 
populated, as for a beam, one finds Pp/pp = —G^ J dp'S{p — p') giving the expected 
depletion of the beam. 5" can of course also account for non-trivial interactions among 
the background particles; it will then have more structure than simply the 5 function. 

Block equations. — An opposite, more "quantum mechanical" limit occurs when 
S depends negligibly on the energy transfer, as can happen, say for a fixed spin, or 
more generally when the energy transfer is small compared to the scale of energies 
associated with processes in the medium. Then ^(A) S'(A) = S'(— A) and wc find 
that the non-linear terms in Eq. (9) disappear, leaving for the matrix structure of 
the collision term 2GpG — GGp — pGG = —[G, [G,p]], which puts the nature of the 
collision term as a kind of double commutator in evidence. The disappearance of the 
non-linear terms is interesting; self-blocking by the ipi docs not reduce the damping. 
Evidently it should not be thought that the particles do not scatter even if they are 
not permitted to change their state in some particular basis. It is in this way, after 
all, that a degenerate Fermi gas exerts a pressure. Blocking effects in the background 
as in Eq. (10), on the other hand, are no^to be ignored. 

Carrying out the integrations in Eq. (9) gives an equation solely in terms of the 
total p (for O, independent of p), leading to a Bloch-hke^ equation. In the two- level 
case this gives for the collision term: 

i J dpdp'S{p - p')npG X (G X P) = -i J dpdp'S{p - p')np|G|2py = -DPt- 

(11) 

T means the direction perpendicular to G, where G — ^{Gq + G ■ r). This is the 
damping term of Eq. (1). The calculations given here verify, by essentially standard 
manipulations, the results of the more intuitive S'-matrix arguments, where D was 
found for two states ("up" and "down") as the real part of (5'up5'^Q^j^ — 1) or in terms 
of phase shifts to lowest order, D ~ (5up — ^down)^? which is indeed ~ {g^p — (?down)^; 
as in Eq. (11). This contains, for example, the old result^ that if one type of particle 
interacts and the other not, that the damping parameter is equal to half of the 
(unblocked) collision rate of the interacting type. 
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Sign of D and Spontaneous Polarization. — D in Eq. (11) has arrived with the 
correct sign — we have damping and not anti-damping. This "arrow of time" arises 
from straightforward calculation — at least for the short time the interaction is on. It 
is a consequence of our starting from an uncorrelated state, which can only become 
more correlated in time. That the damping sign is not entirely trivial is revealed 
by numerical solutions of Eq. (8) where we find when spin and momentum degrees 
of freedom are closely coupled that the relaxation of the momentum can lead to a 
transient polarization for an initially unpolarized situation. 

Neutrino Mixing in a Stochastic Background. — Flavor mixing in a supernova 
would be of significance since it allows other neutrinos to share the large degeneracy. 
In previous work the approach to kinetic and chemical equilibrium was treated in 
the single-particle wavefunction approach^^. Our kinetic equation, with degeneracy 
effects built in, is well suited to this problem. Even though our treatment here 
is not explicitly relativistic, relativistic kinematics can be incorporated in f2p, and 
the absence of the anti-particle degrees of freedom is not critical since these are 
suppressed by the high chemical potential. It is important, however, to include Vp — e 
conversions; this may be done by introducing further bilincars involving the electrons, 
a detailed account will be given elsewhere^^. For applications where the anti-particles 
are important, the formalism should be extended, which raises the interesting question 
of dealing with relativistic fields. 

Understanding "Decoherence" . — We recall^ that in the first quantized formula- 
tion the value of gives the extent to which the eigenstates of G, "up" and "down" 
(two- level case), can interfere. As is reduced the two states "decohere". Similarly 
here, interference terms between operators corresponding to different eigenvalues 
of G are given by the appropriate elements of p and are generally decreased by the 
action of the collision term. Hence we may say that in the collision integral we are 
finding the (rate of) "decoherence". The damping or decoherence rate is hke the 
scattering rate by the background, but with the difference amplitude ~ (g^up — fl'down); 
difference because the scattering must distinguish the two systems in order to "mea- 
sure" and "reduce" them. Hence interactions not distinguishing the states in question 
are not "reducing". The other extreme, that of strong damping, leads through the 
rapid disappearance of Py, to the "freezing" of the other componcnt^'^ and hence 
to the "Turing- Xeno- Watched- Pot" effect^, giving the stabilization^ of the molecular 
isomers. 
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